Spontaneous emission of an atom near a wedge 
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It is a well known fact that the presence of material bodies in the vicinity of an atom affects 
its interaction with the always present quantized electromagnetic field. In this paper, we focus on 
how the spontaneous emission rate of a given excited atom is altered when this atom is placed 
inside a perfectly conducting wedge. We begin by briefly presenting the formalism on which our 
calculations are founded, proceeding then to a long but straightforward calculation of the transition 
rate. We present results for a general atom but, for the sake of simplicity, we narrow them down 
to an effective two-level system in our numerical investigations. The oscillatory pattern for the 
spontaneous emission rate of the atom as we vary its relative position to the wedge as well as the 
phenomenon of suppression are shown in a couple of graphs. 



PACS numbers: 12.20.Ds, 34.20. Cf 

I. INTRODUCTION 

Last year we have celebrated the ninetieth anniver- 
sary of A. Einstein's historical paper on his theory of 
radiationpQ, the final paper of a series that began more 
than a decade earlier with his 1905 work on the quanti- 
zation of the radiation field. By trying to explain how 
atoms could be in thermal equilibrium with the radiation 
field, he introduced the novel concepts of spontaneous 
and stimulated emission of radiation [2], that again out- 
lined how unusual quantum mechanics could be and re- 
heated a dormant discussion about causality in the quan- 
tum scale [3]. However, despite its huge achievements 
this remarkable work had at least one shortcoming, since 
it did not furnish the means for evaluation of the sponta- 
neous emission rate. This had to wait another ten years, 
until P.A.M. Dirac successfully calculated it in a paper 
that most regard as the birth of quantum electrodynam- 
ics [5]. 

Once established the framework of the modern quan- 
tum field theory, a meaningful question to ask is how 
spontaneous emission (SE) was affected by the presence 
of boundaries. That nontrivial boundary conditions could 
modify the SE rate is very clear from the QED stand- 
point, since it is the quantized electromagnetic field that 
causes SE and the field is definitely affected by bound- 
ary conditions This issue was firstly considered by 
E.M. Purcell [7] in the forties, and subsequent experi- 
mental work showed that SE could be strongly modified 
if boundaries were present [8, 9 . The first explicit the- 
oretical results [10l Qj] for parallel conducing plates also 
supported the conclusion that boundary conditions could 



significantly enhance or suppress the SE rate [12] . 

In this work we carry on the analysis on how bound- 
aries may influence the SE process. We investigate the 
behavior of an atom inside a wedge made of conducting 
plates. The reason why we chose such a setup is twofold. 
Firstly, the wedge configuration is a relatively simple sys- 
tem that allows the study of nontrivial curvature effects 
in the zero-point field [I3HI6] . due to the presence of 
the cusp. In addition, the wedge geometry has already 
been used in experiments |17)-although thereabouts the 
wedge aspect was not important, it may probably be in 
the future. Such configuration may also be relevant for 
trapping atoms in excited states, since a suited geometry 
could lead to locations where the potential well disturb- 
ing the atom is deep and the SE rate is small. 

The plan of this paper is arranged as follows. In the 
next section we present the necessary formalism in some 
detail, and show how it leads to a systematic way of eval- 
uating SE rates. In section III we define our problem and 
proceed with the calculations needed in order to bring the 
result to the most convenient form for numerical inves- 
tigations, leaving the fourth section to conclusions and 
perspectives. 

II. TRANSITION RATES OF A SYSTEM 
COUPLED TO A RESERVOIR 

A. General formalism 

Let be a small system S characterized by a time de- 
pendent density matrix operator ps{t) which interacts 
with a large system 1Z which can be treated as a reser- 



2 



voir and characterized by a time independent density ma- 
trix operator p R . The Hamiltonian of the whole system 
S + TZ is H = H s + H R + V, where H S {H R ) is the un- 
perturbed Hamiltonian of S(1Z) with eigenstates |a)(|/z)) 
and with eigenvalues E a (E ll ). By hypothesis, the inter- 
action hamiltonian between S and TZ is 



(i) 



where Sj and Rj are compatible observables associated to 
S and TZ respectively. Assuming that correlations occur- 
ring between S and TZ last a time r c and that variations 
of the observables in S are characterized by a time scale 
Ts, one is able to write a coarse grained rate equation for 
the time evolution of the diagonal elements of ps(t) in 
the base ket \{a}) [HUH], 



dt 



pL (*) 



E 



P S aa(t)T a - 



(2) 



where dt is not infinitesimal but satisfies the inequality 
r c <C dt <C Ts, provided of course that such a relation 
holds for some dt. The factor r c ^ a is defined by 
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5> M ]T|(M, C |!>,a>| 2 . 

S (E^ +E C -E U - E a ) , 



(3) 



where p^ b = (a\ps\b) and is the statistical weight for a 
state \p) of the reservoir. We see that r c _> Q may be inter- 
preted as the probability (per unit time) of a transition 
in the system between the states |c) and \a). Equation 
([3]) is precisely the Fermi's golden rule. 

In order to obtain the exchange energy rates, one must 
consider the time derivative of the average value of the 
hamiltonian of the system. For a given eigenstate \a) of 
the system, a simple application of Eq. |2| shows that 
the net exchange energy rate between S and TZ is 



frequency space [2"U1 |2"T] , 



Cji (w) 



~R 1 \ 

Xjk (w) 



Xjk (w) 



^ R^Rt^ \s {oj + u)n U ) + 6(ui- w M „) J , (5) 

V 

Xjk («) + ixjfc* (w) , (6) 
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where ffl — (p\Rj\v), = (E^ — E v ) /H and V is the 
principal part. The equivalent expressions for the system 
S are 

Cjk M = 7T ^ Sj n S na [ff (u + ui an ) + 6 (w — ui an ) j (9) 

n 



1 



Wan + W 



(10) 



Xj'fc S ' a M = ^ ^ SjnSna [5 (ui an + Uj) - 5 (cj an - a;) j (11) 
n 

with = (a\Sj\b). It is now a matter of straightfor- 
ward algebra to show that Q may be rewritten in terms 
of these new quantities 



d(H. 



S a 



dt 



Q a = Qi r + q: 

I,} V 

2?r 



E 



(12) 

4;rHcgH, (13) 



Q7 = "E 

3,k 
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w*#*(u)C£f M .(14) 



The physical interpretations of Q{ r and Q£ r are simple. 
The former may be understood as the power kicked into 
the system S by the fluctuations of the reservoir (fr), 
while the latter represents the power lost to the reser- 
voir due to fluctuations of the system itself. This last 
term is also called the reservoir reaction contribution (rr) 
since these fluctuations of the system affect the reservoir, 
which then back reacts on the system [TSl HH] . 



= PaaEa = ~ ^ab^a^b , (4) 
b 



where w a (, = (E a — E b )/H. Following the work developed 
by J. Dalibard et al. [THHHI, we now proceed to recast 
the previous result into a more intuitive form that clearly 
exhibits the roles played by S and TZ. Let us begin by 
defining the symmetric correlation function Cj| (ui) and 

the linear susceptibility xfk (^0 °1 ^ ne reservoir in the 



B. Dipole interacting with the radiation field 

Let us now take our system to be a small neutral entity, 
characterized by the set of transition frequencies Wb a , and 
let us assume that it is in contact with the quantized ra- 
diation field, taken as the reservoir. The radiation eletric 
field may be written as 

E (x, t) = ]T (f kA (x) e^4 A + f£ A (x) e— *Va) (15) 

kA 
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where each mode is specified by a wavevector k and polar- 
ization A, with ujk — c|k|. The annihilation and creation 
operators and a^ x satisfy the usual commutation re- 
lations 



t 

k'A 



[okAjfflk'A'] = [ a ikA' a ' 

[<»kA,4'A'] = S] - 



= 0, 

kk'^AA' 



(16) 
(17) 



and the function fkA (x) is the position dependent part of 
the mode kA and carry the information about the bound- 
ary conditions and possible source contributions. 

Let us assume that the dimensions of the system S are 
too small compared to any of its transition wavelengths, 
that is, ao <C 2irc/u!ab, where clq is the largest dimension 
of the system, and also that S is non-relativistic. Then 
the interaction part of the hamiltonian may be approx- 
imated by the dipole interaction of a particle and the 
electric field 



U(x,i) 



: -d • E(x,f) = 



kA 



(18) 



kA j 



where d is the dipole moment operator of the system 
/kA ( x ) 1S t ne ^'-component of fkA (x) and the index 
j = 1,2,3 runs over the directional unitary vectors of a 
given tridimensional coordinate system. Assuming that 
the susceptibility and the correlation function of S are 
diagonal in the appropriate coordinate system, we have 



= -^E 



a ab UJ o.b 
2 



S (uj ab +LU) +5 (uJab - w) ,(19) 



- "S,a i \ 



= a (+) 



E 



a ab Uab 



m 

5 (oj ab -u)—6 {uj ab + uj) 



(20) 



-2|(a|^|6)| 2 /S Wafc . 



where 

Using last equations in ( [T3] ) and (14), making \fi) 
\rikx) in equation ^ (which means a Fock state with n 
photons in the mode kA) and performing the integration 
on uj, we obtain, after some calculations, 



67 = -EE cfea & (*)i^( x 



kA 



Qa = E E C k < W 1/kA w i 2 ( 2 Ka) + 1 

j kA 



(21) 
(22) 



where (n^x) is the average number of photons in the 
mode kA. 



As pointed out before, equation (21) gives the contri 



bution arising from the radiation reaction to the energy 
rates of the system, while equation ( 22 ) gives the corre- 



From ( 21 ) and ( 19 ), it is possible to see that the (rr) con- 



tribution is never positive, which means that it always ac- 
counts for power that is emitted by the system. Since it 
also does not depend on (n^x), we conclude that the (rr) 
contribution comes from spontaneous processes only. On 
the other hand, Equations ( 22 ) and ( 20 ) show that the 



(fr) contribution can have both signs, meaning that Q^J 
could represent either an emitted or an absorbed power. 
In addition, its dependence on (nkA) signals that stimu- 
lated processes also play a role in the fr contribution. For 
the special case where (n^x) = 0, there is no absorbed 
power and the total (spontaneous) radiated power by the 
system is 



Qa = E Co. 



b,j 
b,j 

^7rcY,J2<kabY, k \fL (x)| 2 <5(fc-fc Qb ) , (23) 

j b<a kA 

where k ab = u ab / c and. Let us also note that the emitted 
power at a given permitted transition, Q a b,j, is propor- 
tional to the spontaneous emission rate r^_^ b character- 
istic of this transition. 



III. THE SPONTANEOUS EMISSION OF AN 
ATOM INSIDE A WEDGE 

As far as we know, the method described in the previ- 
ous section was firstly applied in atom-surface problems 
by D. Meschede et al. |22J, where they found the level 
shifts and the radiation rates of an atom near a single 
wall. A few years later, some of the present authors 
revived this method to evaluate thermal corrections to 
the same problem and also to calculate the van der 
Waals interaction between an atom and a perfectly con- 
ducting wedge [24 . In this paper, we intend to continue 
this work by evaluating the spontaneous emission rates 
for the atom-wedge system. 

Let us consider the system described schematically in 
figure 1. It consists of an atom placed inside of a perfectly 
conducting wedge, characterized by a radius R and an 
angle a. The location of the atom is given by its dis- 
tance p to the wedge axis and its declination <p relative 
to the bisector plane, related to the auxiliary angle <fi by 
ip = (j) — a/2. Let us further assume that this atom is 
characterized by a dipole moment d = — er e , where e is 
the fundamental charge and r e is the electron position 
with respect to the center of the atom. It is then clear 
that taking the atom as our system S and the electro- 
magnetic (EM) field constrained by the presence of the 
wedge as the reservoir TZ, we may straightforwardly apply 
the formalism discussed in the previous section. 



According to (21 1- (22), we need the modes of the EM 



sponding contribution coming from the field fluctuations. 



field inside the wedge in order to evaluate the sponta- 
neous emission rate. As they have been already obtained, 
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FIG. 1: The Atom- wedge system. The atom is represented 
by the black circle. 



for instance, in |13j . we will merely quote them 

f ™,n ( x ) = Pqm (lm, n R) (lm,J ~ lk z V t ) • 

[J qm (lm,nP) sin {qm<j)) e - lk > z ] , (24) 

^k,m,n ( x ) = T-kfiqm (Vm,nR) ' 

zxV t [J qm (r}m,np) cos (qmcj)) e~ lk ' z ] , (25) 



where q — w/a, m S N* in (p4|, m G N in ([25]) and 



2 2 r 2 i ; 2 



The quantities 7 m „ and rj m n arc dchncd by 



(26) 



Jqm {lm,nR) = Jqm {Vm,nR} = n € N (27) 



and 



/3gm fa) = -^r- X qm{x) if ™/0. 



/3gm fc) = -jrXqm{x) if TO = . 



(28) 
(29) 



= [J,: 2 (x) + (l-^/x 2 ) J^(x)] (30) 



Vt = pd p + p 



(31) 
(32) 



with K mn0 = 7 m , „ and rc mn i = r? m ,„. Inserting ( 24 1-( 25 ) 



into (23) and setting (nkx) = 0, since we are at zero 
temperature, we get after a little manipulation 



/oo 00 00 
d/Cz X! X K M«A <5( fc af, " J«f m |„ A X -?M ( K |m|nA^) Q#m|, n (p, 0) , (33) 

. . , -oo ___„„ 1 



where A = 0, 1, j — r, </>, z, and 



m— — oo n—1 



Qqm ( K miiOi Pi 

Q g m° («mnO,P, 
Qgm ( K mnOi Pi 
Qgm ( K mnl! ft ^ 

Qqm ( K mnl;P:0 
Qgm ( K mril ) P; 



-7T~Jqln (KmnOP) Sm 2 (gTO0) , 
K mnO 

"~ I 2 J «m ( K mnOP) COS 2 (gm0) , 



(KmnoP) sin 2 (qmcj>) 



= 1 



A 1 



= 1 



K mnlJ K mnlP 



;Jq m (K m 7ilp) sin 2 (qm<f)) 



Jqm («mnlP) COS (gm0) , 



= 0. 



We may now use the generalized Abel-Plana formula [T31 [55] , which is 



oo 1 f°° n 

y^X K mn\R)X qm (K mn \R)f(K mn xR) = - dxf(x) + -Res 

n=l J° 



Jqm (^) 



1 
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e- qmlTl f{e" /2 x) + e qm7rl f(e- l7T/2 x) 



(34) 

(35) 
(36) 
(37) 

(38) 
(39) 



(40) 



where the superscript (A) in the Bessel functions mean their A-th derivative, in order to recast the expression 
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(33) into 



/oo 00 />oo 

dk z ^2 / d V V 3 6 ( k ab - VV 2 + k z) ^' ^' <W + J7 (- R )/°) 

A,,.„ -°° m=-oo"' 



(41) 



where we were able to isolate all the R-dependence of 
Q a into an involved but given function J"(i2, p, 0). As 
we are interested in the wedge without the external cap, 
we must take the R — > oo limit, and here is where the 



rearrangement we made in ( 41 1 shows its quality: it may 
be shown that J-(R — > oo, pTw ~~ ► P3]j leaving us only 
with the first term. 

Although the generalized Abel-Plana formula made 



things simpler, expression (41) still is not very efficient 



for numerical investigations. However, if we rewrite the 
opening angle as a = n/q and restrict ourselves to inte- 
ger values of the parameter q, it turns out that further 
simplification is possible. In this case we can switch Iml 



to m in (41 ) and carry out the m-summation by using an 



addition theorem relating Bessel functions [26] 



oo q—i , , . 

VMW(Me 2 " m * = -]T(-ir /2 e-^+f )Z V ( 2 K psin U+- 

m=-oo * 1=0 ^ \ " 



(42) 



r 



where Z v is a given solution of Bessel equation. As 
slightly different cases of this theorem will apply to each 
polarization, it is convenient to split Q a in its p, <fr and z 
contributions 



(43) 



and consider each one separately. Beginning with the 
z-contribution and substituting (36), (39) and (42) into 
Qa,z, we get 



Qa,. 



H-KC 

^2q 



„ ... 



dk z 

X. Jo 



dyy 3 5(k ab ~ x/y 2 + k 2 z j 



9-1 

E 

Z=0 



J 2y/9sin 



Jo 2ypsm 



<1 



By performing the change of variables y = fcsin#, k z — kcos9 and using the identity |27) 

-7r/2 

J„ (a sin 6) {sm0f +1 (cos8) 2 " +1 d9 = 2T {v + 1) a""" 1 J„ +|i+ i («) , 



(44) 



(45) 



valid when Re [i/] > — l,Re [/x] > —1, we may perform the trivial fc- integration (due to the S- function) to finally put 
( 44 ) in the form 



Q a , z = -~'52to\k ab \Tl+ b i~-G\\ (2\k ab \psm^)-J2 [G\\ (2|fc„ 6 |psm (</» + irl/q)) - G {{ (2\k ab \psm (irl/q))] \ , (46) 

b<a { 1 = 1 ) 



where 



and we have made use of the explicit expressions for 



GJx) 



sin x cos x sin x 



(47) 



() 



J\ii ( x ) an d Jz/2{x)- We have also defined 
r'U b = \MdAb)\ 2 \k ab \\ 



In these cases the m-summations that we need to evaluate 
are 



(48) 



which is nothing but the spontaneous emission rate be- 
tween levels a and b for the j-polarization in free space, 
assuming of course that a > b. 

The calculations for Q a . p and Q a ,4> are a little bit more 
involved but rather analogous to the one we just showed. 



m2j qm( K P) 



cos 2 (qmcj)) 
sin 2 (qmcj)) j ' 

cos 2 (qm(f>) 



sin (qmcj)) 



(49) 
(50) 



The first ones, shown in (49), may be carried out by considering the particular case of (421 in which Z u+qm (np) 
Jqm(Kp) and differentiating it twice with respect to the variable cj>. After some algebraic manipulations we get 



q-i r 



q 2 1 

5^ m 2 J 2 qm (up) cos {2qm(f) = - ^ 

q 1=0 



K 2 p 2 



m— — oo 



Jo (2np sin ipi) cos ij>i 



2 Ji (2fcp sin ■;/>;) 



2np sin tpi 



(51) 



where ipi = cj> + trl/q, which immediately leads to 

9-1 



m 2 J qm (up) sin 2 {qmcj)) = — [—Jo (2fcpsin^) cos 2 V>z + Jo (2«p sin 7tZ/q) cos 2 nl/q 
' Lq 



m= — oo 



1=0 



Ji {2nps\aipi) Ji (2Kpsin7r//g) 



2 up sin ^i; 



2«;p sin irl/q 



(52) 



9-1 



— ' 1 

m 2 J qm (up) cos 2 {qmcj)) = — [Jo (2«;psm ?/>;) cos 2 "0; + Jo (2Kpsm.nI/q) cos 2 7r//g 



Ji (2Kpsin^;) Ji (2Kpsimrl/q) 
2Kpsin^; 2Kpsmirl/q 
The remaining summations in (1501) may be evaluated by making use of the identity 



(53) 



Kp 



Jqm ( K P) Jqm («p) ~ 9 Kp [J qm («p) Jj m+ 1 (ftp)] 



t2 p + 1 T , . T , . 

+ Jqm i K P) Jqm («W Jj m+ 1 (ftp) , 



(54) 



plus the particular cases of (42 1 for Z v + qm {np) = J qm (np), J qm +\{Kp) and its derivatives with respect to variables p 
and (p. This yields 



9-1 



J qm (kp) sin 2 (qmcj)) — — [Jo (2«p sin ipi ) sin 2 ^ — Jo (2Kpsmirl/q) sin 2 nl/q 

„— 9 



i=o 



9-1 



Ji (2Kpsin-0;) Ji (2Kpsin7r//g) 



2np sin -0; 



2/ipsin 7rZ/g 



(55) 



J 9 ' ? 2 (rep) cos 2 (qmcj)) = — [—Jo (2repsin^>z) sin 2 — J (2«;psm7rZ*/g) sin 2 ttZ/q 



m— — oo 



Ji (2«;psin'0;) Ji (2np sin nl/ql) 



2np sin ^ 



2/tp sin7rZ/g 



(56) 



By inserting equations (52) - (56) into the expressions for Q a , P , Q a ,rj> and following the same procedure that led to 
( 46 1 , we finally get 



9-1 



2a,0 = ~^c|fc a6 |r o O |-^(2|* a 6|p,0)-^[^(2|fe a 6|p,^ + 7rZ/g)+^(2|fc a6 |p ) 7rZ/ g )] L (57) 

b<a { 1=1 ) 
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where 



H^(x,if)) = Gy (xsin^) sin 2 if) 

+2Gj_ (x sin if)) cos 2 V (59) 
H p (x,ip) = Gy (xsin?/>) cos 2 if> 

+2Gj_ (x sin t/>) sin 2 if) , (60) 



As we now have well suited expressions for the decay- 
ing rates, we may proceed with some numerical investi- 
gations. In order to simplify our discussion, let us as- 
sume that our atom has only one dominant transition, so 
that we can neglect all terms in the summation of states 
present in ( 46 1-( 58 ) except for this dominant term. 



and also 



G x (x) 



cos x sin x 



(61) 




FIG. 2: The spontaneous emission for the z-polarization in a wedge of 7r/3 radians. The axes X, y and Z follow the same 
conventions of figure 1. The plates are not displayed but the plot itself is very suggestive - they would be sharply cutting the 
two highest slopes in half along the z direction. Their aspect can be easily inferred from the diffuse V-shape drawn by those 
slopes. 

I 



The first thing that we would like to point out is how 
differently each polarization responds to the presence of 
the wedge. By one hand, as one can see from Fig. [2] 
the contribution of the z-polarization for the sponta- 
neous emission inside a wedge is greatly suppressed in the 
vicinity of the plates. However, by taking a look at the 



equivalent graph for the ^-polarization shown in Fig. [3j 
one concludes that just the opposite occurs, which means 
that this contribution is enhanced near the plates. This 
effect may seem curious initially, but it is readily under- 
stood by recalling the boundary conditions imposed by 
the wedge on the electric field, given by ^ x E| s = 0. 




This implies that very close to one of the plates and rel- 
atively away from the cusp the electric field is approxi- 
mately perpendicular to the mirror surface. Since we are 
working in the dipole approximation, we conclude that 
the parallel component of the dipole no longer couples 
to the field while the perpendicular component coupling 
is maximized. This explains both the vanishing of Q a , z 
(and also of Q a ,p) and the enhancement of Q a ^ near the 
plates. We should mention that the previous reasoning 
is not valid when that atom is too close to the cusp, and 
we see accordingly that in such a case the SE vanishes 
for both polarizations. This is, of course, a direct conse- 
quence of the vanishing of the electric field at the cusp, 
as it may be seen from expressions ( 24 1 and ( 25 1 . 



Intuitively, we should expect that a wedge character- 
ized by very small angles should mimic the behavior of 
two parallel plates, provided the atom is not too close to 
the cusp. The fact that we have checked numerically that 
our expressions coincide with the known results for par- 
allel plates in the limit a — > 0, p — > oo , pa const., en- 
couraged us to consider less extreme situations, in which 
a is small but not infinitesimal. An interesting thing 
happens when the setup under consideration is like the 
one depicted in Fig. 4, in which an atom with Q a ,4> = 
is placed in the X axis at a given distance from the cusp 
p = x. If a <C 1 it is reasonable to approximate the wedge 
by two parallel plates, but the distance between these ef- 
fective plates changes as we move the atom along the X 
axis, like illustrated in Fig 5. So, by considering an atom 
in different positions at the bisecting plane of an acute 
wedge, we are approximately also looking to the situation 
in which the atom is halfway between two parallel plates 
and these plates are put at different separations. In Fig. 



move the atom along the X axis for different values of q. 
We see that all curves are practically zero up to approxi- 
mately k a bX = q, where they rise abruptly and then start 
oscillating around a decaying mean value. This is more 
clearly understood in terms of effective parallel mirrors, 
as we shall see. 




6 we show the behavior of Q a 



Qa,z + Qa,n &S we 



Y 



FIG. 4: An atom localized in the bisecting plane of a wedge 



The perpendicular component of the wave vector char- 
acterizing the modes between two plates is discrete and 
given by mr/5, where S is the distance between the plates. 
When k a fjX < q, we have necessarily S < A a b/2 (remind- 
ing that S ~ xa = xir/q and X a b = 2ir/k a b). This im- 
plies that all the modes with n > are more energetic 
then the a — > b transition, making their excitation by 
the dipole impossible. That said, we conclude that the 
only way in which the dipole could emit would be to 
populate the mode n = 0, but that is also impossible 
since the zero mode is polarized perpendicularly to the 
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plates (T2], not coupling to the dipole at all [28 . The 
situation changes dramatically once k ab x ~ q, since the 
availability of at least one mixed (with perpendicular and 
parallel components) mode with E moc i e < huiab allows the 
atom to emit, explaining the abrupt jump at that point. 
Based on the same reasoning, we should expect similar 
jumps each time the distance between the atom and the 
cusp reaches a multiple of X a b/2 (or equivalently, if k ab x 
reaches a multiple of q), but as we see from the curves 
corresponding to q = 60 and q = 90, there are jumps at 
k ab x ~ 3g but not at k ab x « 2q. This may seem strange 
at first sight, but it is a direct consequence of the struc- 
ture of the allowed modes between two mirrors. As we 
said before, these modes are characterized by an integer n 
present in the perpendicular component of the respective 
wave vectors. By using the explicit form of these modes 
[12], it is possible to show that the ones associated with 
even integers have a node at 5/2, producing a vanishing 
electric field there. Since our atom is effectively located 
at 6/2, we conclude that it feels only the influence of the 
odd modes, and that is the reason for the jumps only at 
k a bX equal to odd values of q. 




I 




FIG. 5: Effective mirrors at two different positions of the 
atom. We note that the distance between the effective plates 
becomes larger as the atom moves away from the cusp. 



We shall close this work by considering a limiting case 
of our expressions, in which we put a = n and the wedge 
degenerates into a single plate. Then q = 1 and by eval- 



uating (46), (57) and (58) at = tt/2 we get 

QaA\ = Qa.x + Qa,z = -\hk ab c{T X a ^ b + I 



a^b) 



2 _ sin(2|fc ab |x) 

3 2\k ab \x 



cos(2|fc ab |a:) 
(2\k ab \xY 



Qa,y 



-hk ab cT y a ^ b 



sin(2|fc a6 |x) 
(2\k ab \x) 3 

_ cos(2|fc ab |a;) 

{2\k ab \xf 
sin(2|fc ab |x) 
(2|fc a6 |x) 3 



,(62) 



(63) 



which is in precise agreement with the known result from 
the literature [HE]. 

IV. CONCLUSION AND PERSPECTIVES 

In this paper we investigated how the spontaneous 
emission rate of an atom is affected by the presence of 
a conducting wedge. We used a very general method, 
which allows for the evaluation of van der Waals and ex- 
cited potentials as well as spontaneous emission rates. 
Firstly we applied it to the problem of an atom inside a 
mirror wedge. In the following section we showed that, 
although for a wedge of an arbitrary angle a the expres- 
sions are considerably involved, when q € N a lot of sim- 
plification was possible and we arrived at manageable ex- 
pressions for SE rates. We then obtained the contribution 
of each polarization and plotted two of them separately. 
Note that the behavior of the SE rate in a wedge shows 
the usual oscillatory pattern. 

As expected, the phenomenon of suppression of the SE 
rate can also occur for an atom inside a wedge. As we 
showed graphically, for very small angles of the wedge, 
the plates forming it behave effectively as they were two 
parallel plates. Hence, for polarizations parallel to these 
effective plates, namely, polarizations p and z, there are 
configurations for which the atom will not decay at all. 

There are several possible generalizations for this work, 
the most obvious one being how to bring some real effects 
into consideration, like finite conductivity and tempera- 
ture corrections, which fall out of the scope of this paper 
but serve as a guide for future research. 
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